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Abstract 

The DFT/vdW-WF2sl method, recently developed to include the van der Waals interactions 
in the Density Functional Theory and describe adsorption processes on metal surfaces by taking 
metal-screening effects into account, is applied to the case of the interaction of Xe and graphene 
with a transition-metal surface, namely Ni(lll). In general the adsorption of rare-gas atoms on 
metal surfaces is important because is prototypical for physisorption processes. Moreover, the 
interaction of graphene with Ni(lll) is of particular interest for practical applications (efficient 
and large-scale production of high-quality graphene) and, from a theoretical point of view, is 
particularly challenging, since it can be described by a delicate interplay between chemisorption 
and physisorption processes. The first-principles simulation of transition metals require particular 
care also because they can be viewed as intermediate systems between simple metals and insulating 
crystals. Even in these cases the method performs well as demonstrated by comparing our results 
with available experimental data and other theoretical investigations. We confirm that the rare 
gas Xe atom is preferentially adsorbed on the top-site configuration on the Ni(lll) surface too. 
Our approach based on the use of the Maximally Localized Wannier Functions also allow us to 
well characterize the bonds between graphene and Ni(lll). 
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I. INTRODUCTION 


The interaction of graphene with transition metal surfaces is of great importance because 
the growth of graphene on these substrates is probably at present the most perspective way 
for the large-scale preparation and production of high-quality graphene.-i^ In particular, 
the Ni(lll) surface is very interesting, being the closest matched interface with respect 
to graphene of all transition metals (with a lattice mismatch of about 1%), and has been 
therefore intensively investigated.-^— In fact, the close lattice match allows graphene to 
adapt itself to the Ni(lll) lattice and enables the formation of a well-ordered p(l x over¬ 
structure on the Ni(lll) surface, which is much simpler than the complex Moire patterns 
commonly found in other transition metal surfaces.-^- Moreover, given the ferromagnetic 
properties of Ni, graphene on Ni surface has also been proposed as a promising spin-filtering 
device needed in spintronics.— Due to the strong interaction between graphene tt and Ni 3d 
electrons, the graphene electronic structure is heavily modified if compared to the electronic 
structure of free-standing graphene. The most important change is the gap opening resulting 
from sublattice symmetry breaking. 

From a fundamental point of view, the nature of the bonds at the graphene-metal inter¬ 
face has not been completely elucidated yet. Basically, by taking the C-C covalent bonds as 
a reference, the interaction of graphene with a metal substrate can be classified as ’’weak” 
(with Ir(lll), Au(lll), Cu(lll), Pt(lll), Al(im and Pd(lll)) or ’’strong” (with Ni(lll), 
Co(OOOl), Rh(lll), and Ru(OOOl)), see refs. ^,21 and references therein. Clearly this be¬ 
havior depends on the atomic and electronic structures of the different substrates which 
can substantially modify the properties of the adsorbed graphene with respect to those of a 
free-standing graphene layer. In particular, transition metals can be viewed as intermediate 
systems between simple metals and insulating crystals because their conduction electrons 
are of both free-electron (s- and p-like) character and of localized (d-like) character, so that 
they are neither totally free nor completely bound. Moreover, the d-like electrons of the hlled 
bands are not rigidly bound to the ionic cores and therefore must be somehow included in 
the description of screening effects. 

From a theoretical point of view the interaction of graphene with Ni(lll) represents a 
particularly challenging system, since, both an accurate description of the metallic surface 
and of the nonlocal correlation effects is needed. Moreover, in the most favored adsorption 
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configurations, van der Waals (vdW) effects are very important although the graphene-metal 
interaction is relatively strong and the equilibrium distance of about 2.1 A is much shorter 
than typical distances observed between fragments bonded by pure vdW interactions, thus 
suggesting that a delicate interplay between chemisorption and physisorption exists.—*^ As a 
result, the computed adsorption energies and equilibrium distances, even by adopting vdW- 
corrected Density Functional Theory (DFT) approaches, exhibit a disappointing wide spread 
of values.— For instance, using the Generalized Gradient Approximation (GGA) functionals 
one has a very weak binding or no binding at all (see ref. Il9| and references therein), while 
the Local Density Approximation (LDA) calculations!^ lead to a pronounced minimum with 
an adsorption energy of about -200 meV per C atom, at an equilibrium distance of 2.00 A a 
little shorter than the experimental value of 2.11±0.07 A.- Actually, the deep LDA minimum 
is a true artifact, being related to the wrong decay of the Kohn-Sham potential and orbitals 
from the surface into the vacuum. In any case the relatively weak bonding confirms that 
the adsorption should not be viewed as a traditional (covalent) chemisorption characterized 
by typical adsorption energies of 0.5-2 eV, although the short graphene-Ni distances clearly 
indicate a covalent interaction already at the level of semilocal GGA functionals. 

The adsorption of graphene on Ni(lll) has been also investigated on the basis of the 
adiabatic-connection fluctuation-dissipation theorem in the random phase approximation 
(RPA).— At the moment this probably represents the more accurate calculation, which 
can therefore be taken as a reasonable reference database, being particularly suited to de¬ 
scribe intricate bonds with mixed covalent and dispersive character,—*^ although one should 
remember that RPA typically exhibits a tendency to underbind.— In spite of a significant 
hybridization found between the graphene vr orbitals and Ni <7^2 states at a binding dis¬ 
tance of about 2.2 A, the computed RPA adsorption energy is still in the range of a typical 
physisorption interaction (about -70 meV per C atom). An important contribution to the 
energy is related to a decrease in the exchange energy resulting from the adsorption-induced 
lower symmetry in the graphene layer. Interestingly, the RPA calculations do not only pre¬ 
dict one single minimum but also a second minimum at a distance of 3.3 A. This second 
minimum, with an adsorption energy of -60 meV, that is only slightly less bound than the 
most favored configuration, is found at a distance typical for vdW-adsorbed graphene on 
transition metals.— 

This complex scenario has been rationalized as follows:— at large distances, namely up to 
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2.8 A, the electronic band structure of graphene on Ni(lll) is hardly modihed with respect 
to that of the free-standing, the exchange interaction is purely repulsive, and the correlation 
follows essentially a vdW-like behavior; instead, at distances shorter than 2.8 A, the graphene 
band structure starts to be modihed with a hybridization setting in, and, most importantly, 
with a breakup of the symmetry between the fop-site and hollow-site C atoms. As a delicate 
balance between exchange and correlation energy the total energy is smooth, with a slight 
barrier between the physisorption minimum, characterized by a graphene band structure 
that is hardly modihed compared to free-standing graphene, and another ’’chemisorption” 
minimum, where the graphene band structure is strongly modihed compared to the free¬ 
standing layer. The physisorption minimum originates entirely from vdW forces, whereas 
the chemisorbed minimum can be attributed to both orbital hybridization and week vdW 
forces. The slight energy barrier between physical adsorption and chemisorption state implies 
that the graphene can be decoupled easily from Ni(lll) surface, which is also evidenced in 
chemical vapor deposition (CVD) growth of graphene on Ni surface.— Besides vdW ehects 
at large distances it is therefore also important to describe exchange interaction at a short 
distance where the correct description of Pauli repulsion and hybridization is crucial.— 

Since graphene has two carbon atoms in the unit cell, while the Ni(lll) surface has only 
one, six different structural models are possible, namely top-fcc, top-hep, bridge-top, fcc-hcp, 
bridge-fee, and bridge-hep. In the top-fee and top-hep models, one C atom sits on top of 
the Ni atom while the other C atom occupies, respectively, the fee and hep hollow sites. 
The bridge-top conhguration has the carbon-carbon bond of graphene on top of a surface Ni 
atom. In the fee-hep model, both C atoms are in the 3-fold hollow sites (see Fig. 1). 

Recently, the structural properties of graphene/Ni(lll) have been experimentally 
investigated^i by a combination of low-energy electron diffraction (FEED), x-ray photo¬ 
electron spectroscopy (XPS), and angle-scanned photoelectron diffraction (PED). XPS data 
indicate that graphene interacts strongly with the topmost Ni layer and diffraction data 
show that graphene is deposited commensurably with the underlying Ni surface atoms. 
Interestingly, related hrst-principles calculations based on DFT show^i that the total ener¬ 
gies of the top-fee and bridge-top structures are nearly degenerate, which corroborates the 
observecy^i^ coexistence of those phases. A depletion of charge density in Ni-d 5,2 of topmost 
Ni atoms is also found and an increase of that in p^ of on-top C atom in the top-fee and 
top-hep structures, which indicates that the electrons transfer from substrate to adsorbed 
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graphene.— 

From a theoretical point of view, it is well accepted that the widely-used LDA and semi¬ 
local GGA functionals fail to capture vdW forces, which are important in rare-gas/metal 
and graphene/metal systems. As mentioned above, the most reliable RPA calculations22iM 
show that the binding in graphene/Ni system is a delicate balance between covalent and 
dispersive interactions. However, the RPA calculation is computationally expensive (scaling 
as with system size) and thus hardly accessible for large systems. The cheaper DFT-D 
approach^^ with semi-empirical corrections for vdW interactions seems to typically give a 
reasonable prediction of the adsorption properties. However, its empirical character->2i^— 
casts doubts on the full reliability and accuracy in many applications, particularly in metal 
system where the atom-based description, implied in the use of the Cq coefficients adopted 
in DFT-D, is rather questionable. 

In the last few years several practical methods have been proposed to make DFT cal¬ 


culations able to accurately describe vc 
recent review, see, for instance, refs. l37 


W effects at a reasonable computational cost (for a 


39l ). We have developed a family of such methods, 
all based on the generation of the Maximally Localized Wannier Functions (MLWFs),— suc¬ 
cessfully applied to a variety of systems, including small molecules, water clusters, graphite 
and graphene, water layers interacting with graphite, interfacial water on semiconducting 
substrates, hydrogenated carbon nanotubes, molecular solids, and the interaction of rare 
gases and small molecules with metal surfaces.—^— In particular, the DFT/vdW-WF2sl 
method, presented in ref. ImI, has been specihcally developed to take metal-screening effects 
into account and has been applied to the study of the adsorption of rare gases and small 
molecules on different metal surfaces, namely Al(lll), Gu(lll), and Pb(lll). 

Here we apply the DFT/vdW-WF2sl approach to investigate the interaction of the Xe 
atom and of graphene with the Ni(lll) surface. Our results will be compared to the best 
available, reference experimental and theoretical values, and to those obtained by other DFT 
vdW-corrected schemes, including PBE-|-D^ (a particular DFT-I-D scheme), vdW-DF,— 
vdW-DF2,— rVVlO,— and by the simpler Local Density Approximation (LDA) and semilocal 
GGA (in the PBE flavor—) approaches. In the PBE-I-D scheme DFT calculations at the 
PBE level are corrected by adding empirical Cq/R^ potentials with parameters determined 
by fitting accurate energies for a large molecular database, while in other methods, such as 
vdW-DF, vdW-DF2, and rVVlO, vdW effects are included by introducing DFT nonlocal 
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correlation functionals. 


II. METHOD 


Here we briefly review the DFT/vdW-WF2sl method; additional details can be found 


in refs. 


50 


51 


Basically, the scheme relies on the well known London’s expressiort^ where 


two interacting atoms, A and B, are approximated by coupled harmonic oscillators and the 
vdW energy is taken to be the change of the zero-point energy of the coupled oscillations 
as the atoms approach; if only a single excitation frequency is associated to each atom, 
ujb, then 


rpLondon _ __ ^ \ 

2m‘^ UAUJBioJA +OJb) Rab 

where Za,b is the total charge of A and B, and Rab is the distance between the two 
atoms (e and m are the electronic charge and mass). 

Now, adopting a simple classical theory of the atomic polarizability, the polarizability of 
an electronic shell of charge eZi and mass mZi, tied to a heavy undeformable ion can be 
written as 


a,: ~ 




mu, 


2 • 


( 2 ) 


Then, given the direct relation between polarizability and atomic volume,— we assume 
that Q!j = ■ySf, where 7 is a proportionality constant, so that the atomic volume is expressed 
in terms of the MLWF spread. Si. Rewriting eq. ([T]) in terms of the quantities dehned 
above, one obtains an explicit expression for the Cq vdW coefficient: 


“ 2 (x/zisf + ' 

The constant 7 can then be set up by imposing that the exact value for the H atom 
polarizability [an =4.5 a.u.) is obtained (of course, in the H case, one knows the exact 
analytical spread. Si = Sh = \/3 a.u.). 

In order to achieve a better accuracy, one must properly deal with intrafragment MLWF 
overlap (we refer here to charge overlap, not to be confused with wave functions overlap). 
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This overlap affects the effective orbital volume, the polarizability, and the excitation fre¬ 
quency (see eq. ([ 2 ])), thus leading to a quantitative effect on the value of the coefficient. 
We take into account the effective change in volume due to intrafragment MLWF overlap 


by introducing a suitable reduction factor ^ obtained by interpo 
cases of fully overlapping and non-overlapping MLWFs (see ref. 
the following expression for the Cq coefficient: 


ating between the limiting 


501) . We therefore arrive at 


where represents the ratio between the effective and the free volume associated to 
the A-th and B-th MLWF. 

Finally, the vdW interaction energy is computed as: 


i<j o' 

where f{Rij) is a short-range damping function, which is introduced not only to avoid the 
unphysical divergence of the vdW correction at small fragment separations, but also to 
eliminate double countings of correlation effects (in fact standard DFT approaches are able 
to describe short-range correlations); it is defined as: 


The parameter Rg represents the sum of the vdW radii Rg 
adopting the same criterion chosen above for the 7 parameter) 


( 6 ) 

ifVdW ^ 


rfvdW _ T^vdW 

where is the literature^ ( 1.20 A) vdW radius of the H atom, and, following Grimme et 
al.,— a ~ 20 (the results are almost independent on the particular value of this parameter). 
Although this damping function introduces a certain degree of empiricism in the method, we 
stress that a is the only ad-hoc parameter present in our approach, while all the others are 
only determined by the basic information given by the MLWFs, namely from first principles 
calculations. 

To get an appropriate inclusion of metal screening effects a proper reduction coefficient 
is included by multiplying the Cq /R^j contribution in eq. ([5]) by a Thomas-Fermi factor: 
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Jtf = where is the Thomas-Fermi screening length relative to the electronic 

density of an effective uniform electron gas (’’jellium model”) describing the substrate (see 
discussion below), Zg is the average vertical position of the topmost metal atoms, and zi is 
the vertical coordinate of the Wannier Function Center (WFC) belonging to the substrate 
(/ = i if it is the i-th WFC which belongs to the substrate, otherwise I = j); the above Jtf 
function is only applied if zi < Zg, otherwise it is assumed that Jtf = 1 (no screening effect). 

A. Computational details 

We here apply the DFT/vdW-WF2sl method to the case of adsorption of Xe and 
graphene on the Ni(lll) surface. All calculations have been performed with the Quantum- 
ESPRESSO ab initio package^ (MLWFs have been generated as a post-processing calcula¬ 
tion using the WanT package^). Similarly to our previous studies^*^ we modeled the metal 
surface using a periodically-repeated hexagonal supercell, with a (\/3 x \/3)R30° structure 
and a surface slab made of 15 Ni atoms distributed over 5 layers (repeated slabs were sep¬ 
arated along the direction orthogonal to the surface by a vacuum region of about 25 A to 
avoid signihcant spurious interactions due to periodic replicas), considering the experimental 
Ni(lll) lattice constants. The Brillouin Zone has been sampled using a 6 x 6 x 1 fc-point 
mesh. Above one of the surface Ni layers we add a Xe atom or a single graphene layer. We 
remind that the Ni(lll) surface and the graphene lattice constants are, respectively, 2.489 
and 2.46 A, which corresponds to a lattice misht of about 1.2%, so that one can assume 
that the graphene layer lattice vectors are commensurate with those of the Ni(lll) surface. 
In this model system, with a Xe atom per supercell, the coverage is 1/3, i.e. one adsorbed 
adparticle for each 3 metal atoms in the topmost surface layer. The (\/3 x \/3)i?30° struc¬ 
ture has been indeed observed^ at low temperature by FEED for the case of Xe adsorption 
on Cu(lll) and Pd(lll) (actually, this is the simplest commensurate structure for rare-gas 
monolayers on close-packed metal surfaces and the only one for which good experimental 
data exist), and it was adopted in most of the previous ab initio studies.—^— The metal sur¬ 
face atoms were kept frozen (of course after a preliminary relaxation of the outermost layers 
of the clean metal surfaces) and only the vertical coordinate (perpendicular to the surface) 
of the Xe atom or the graphene layer was optimized, this procedure being justihed by the 
fact that only minor metal surface atom displacements are observed,— and relaxation 


effects are estimated to be small.— Also in the case of graphene, which seems to chemisorb 
on Ni(lll), the buckling between nonequivalent carbon atoms is negligible, indicating that 
the adsorbed graphene layer is quite smooth.— Given the ferromagnetic character of Ni, spin 
polarization was taken into account. 

We have carried out calculations for various separations of the Xe atom and the graphene 
layer, with Xe adsorbed on the top and hollow (on the center of the triangle formed by the 3 
surface metal atoms contained in the supercell) high-symmetry sites, and considering the top- 
fcc, bridge-top, and fcc-hcp conhgurations for graphene on the metal substrates (see above). 
For a better accuracy, as done in previous applications on adsorption processes,— 
we have also included the interactions of the MLWFs of the adsorbate not only with the 
MLWFs of the underlying surface, within the reference supercell, but also with a suffi¬ 
cient number of periodically-repeated surface MLWFs (in any case, given the R~^ decay 
of the vdW interactions, the convergence with the number of repeated images is rapidly 
achieved). Electron-ion interactions were described using norm-conserving pseudopotentials 
by explicitly including 10 valence electrons per Ni atom. We chose the PBE^ reference DFT 
functional, which is probably the most popular GGA functional. 

The choice of a suitable effective value for the uniform electron density n (or, equiva¬ 
lently, the dimensionless ’Vs” parameter, = (S/dvm)^'^^/ oq, where oq is the Bohr radius) 
associated to the metal substrate to get the Thomas-Fermi length dehned above to 
describe screening effects, deserves a specihc comment since we are considering transition 
metals where, differently form simple metals, the electrons cannot be assumed to be totally 
free. To address this issue we used the list of values for the effective free-electron density 
parameter provided by Perrot and Rasolt^ for transition metals. These values were ob¬ 
tained by considering the problem of appropriately describing the ground state properties 
(particularly those of arbitrary defects) of the mobile part of the electron fluid in the tran¬ 
sition metals and dehning the response of the ’’free” mobile part of the electron fluid using 
an effective of a uniform electron gas. The recipe is based on the concept of the metallic 
response to external low-symmetry perturbations and the obtained picture is that, as ex¬ 
pected, the s and p electrons carry the main response to an external potential while the d 
electrons remain largely unpolarized. From the parameter can be evaluated using the 
following formulas:— 
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ttf — ^/krF — ’ (^ 9 ) 
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where and kxF are the Fermi and Thomas-Fermi wavevector, respectively, and g{eF) = 
rnkF/h^F"^ is the density of levels at the Fermi energy eF in the jellium model. Using the 
Tg = 2.14 value suggested for Ni by Perrot and Rasolt,— one obtains a Thomas-Fermi 
screening length = 0.496 A. 


III. RESULTS AND DISCUSSION 

In Table I we report the binding energy and the equilibrium distance for Xe on Ni(lll). 
Note that all the methods predict that the top configuration is favored with respect to the 
hollow one, in line with the observed general tendency of Xe and Kr for adsorption on 
metallic surfaces in the low-coordination top sites.— This behavior is attributed^i^ 
to the delocalization of charge density that increases the Pauli repulsion effect at the hollow 
sites relative to the top site and lifts the potential well upwards both in energy and height. 

The experimentally measured adsorption energy, Ea, often includes not only the interac¬ 
tion of adparticles with the substrate but also lateral, vdW interfragment interactions.—*^ 
Therefore sometimes it is more appropriate to compare experimental data with the quantity 
Ea which can be related to E^, by:— 

Ea = Eh (El — Ef) , (11) 

where Ei is the total energy (per particle) of the 2D lattice formed by the adparticles only 
(that is as in the adsorption configurations but without the substrate and including vdW 
interfragment corrections when vdW-corrected methods are used) and Ef \s the energy of 
an isolated (free) adparticle. 

Experimental estimates of the adsorption energy for Xe on Ni(lll) range from -370 to -180 
meV.— These data were obtained from measurements of the isosteric heat of adsorption^ 
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and optical differential reflectance stndies of adsorption and desorption of Xe on Ni(lll).— 
In the latter investigations the overlayer structure at the saturation coverage (1 Xe atom 
per 3 top Ni atoms) appears commensurate with the underlying Ni(lll) lattice, which is 
consistent with the fact that the Xe-Xe separation on Ni(lll) is 4.31 A, which is compressed 
from the bulk Xe equilibrium value of 4.34 A by only 0.7%. 

Our computed binding energy of Xe on Ni(lll) is similar to that obtained by the rVVlO 
functional, while it is smaller (in absolute value) than the values predicted by LDA and 
PBE+D and larger than the vdW-DF2 value (and obviously the PBE one since this func¬ 
tional does not properly include vdW interactions). Interestingly, by taking the Xe-Xe 
interaction into account to estimate Ea, our DFT/vdW-WF2sl value is compatible with the 
experimental estimates, while instead LDA and PBE-I-D appear to signihcantly overbind. 
To make the comparison with the experiment more accurate one must point out that, al¬ 
though Xe on Ni(lll) is considered a typical case of physisorption process, a net transfer of 
electronic charge is observed from the Xe atom to the metal, which leads to a surface-induced 
weakening of the Xe-Xe interaction resulting from the Coulomb repulsion due to the charg¬ 
ing of the adatoms.— As a result, the absolute value of the correction factor | {Ei — Ef) \ = 88 


meV in eq. flTT]) is probably overestimated; in fact, according to ref. 


82 


, the difference be¬ 


tween Ea and E}, should be in the range of 30-40 meV. Considering this further correction, 
our DFT/vdW-WF2sl estimate for Ea would be about -320 eV, always in line with the 
experimental values. 

In Table II we report our computed binding energies and equilibrium distances for 
graphene on Ni(lll), compared to data obtained by other methods and literature experi¬ 
mental and theoretical estimates (in this case much more data are available than for Xe on 
Ni(lll)). The binding energy curves for graphene on Ni(lll), obtained using our DFT/vdW- 
WF2sl method and relative to the three considered conhgurations, are plotted in Fig. 2. 

By focusing on the top-fcc conhguration, as discussed above, although not necessarily 
exact,— the RPA estimate^'^'^ can be assumed to be the ’’best” theoretical estimate to be 
taken as a meaningful reference. As expected, the PBE functional, which does not include 
genuine vdW effects, dramatically underestimates the binding so that the binding energy 
is positive which means that the system is actually unbound, although the equilibrium 
distance, corresponding to a local minimum in the binding energy, agrees well with the 
reference estimates. However, even the performances of some vdW-corrected functionals 
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are not satisfactory by stressing the difficulty of properly describing this system. In fact, 
for instance, rVVlO largely overestimate the equilibrium distance and, by comparison with 
the RPA curve, it appears to completely miss the first (chemisorption) minimum. The 
vdW-DF and vdW-DF2 methods underestimate the binding energy but (above all) predict 
a much too large equilibrium distance. Other approaches (see Table II) give reasonable 
equilibrium distances, although the estimated binding energies are significantly scattered, 
with the optB86b2^ and optB88vdW— methods which appear to give the best performances. 

For a better assessment of the different approaches, one should remember that RPA tends 
to underbind:— for instance, in the case of the interaction of CO with the Cu(lll) surface, 
RPA remarkably predicts the correct favored adsorption site, however it underestimates^ the 
adsorption energy by 70 meV, corresponding to a significant error of about 15% with respect 
to the both the experimental reference value and high-level quantum chemistry (CASPT2) 
calculations. By taking this observation into account and looking at both the binding en¬ 
ergy and the equilibrium distance, the performances of DFT/vdW-WF2sl are rather good, 
being also better that those of the semiempirical DFT-D2 and DFT-D3 approaches: in 
fact, DFT-D2 tends to overbind, while DFT-D3 evidently overestimates the equilibrium 
distance. Interestingly, the DFT/vdW-WF2sl results are similar to those obtained by the 
PBEsol functional, but for the fcc-hcp configuration. This behavior is understandable since 
PBEsol is a semilocal GGA functional^ (a revision of PBE to better describe solid state 
and surface systems) which does not properly include vdW interactions and therefore cannot 
reproduce the energy minimum of the fcc-hcp structure that is essentially determined by long- 
range vdW effects. Our DFT/vdW-WF2sl results indicate that the top-fee configuration 
is only marginally more favored than the bridge-top^ being in line with recent experimental 
measurements,— >2i which suggest that the total energies of the top-fee and bridge-top struc¬ 
tures are nearly degenerate since the two phases are observed to coexist. Note that LDA and 
PBEsol, which do not take vdW effects into account, differently from the vdW-corrected 
DFT/vdW-WF2sl, DFT-D2, and rVVlO methods, indicate the bridge-top as energetically 
more favored than top-fcc] although this result should not be overemphasized, given the rel¬ 
atively small energy differences between the two configurations, nonetheless it suggests that 
a proper inclusion of vdW interactions can be of importance even when covalent interactions 
are involved, as already observed elsewhere.— Interestingly, the binding-energy curve of the 
bridge-top configuration in Fig. 2 exhibit a shallow minimum at about 3.2 A in addition 
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to the first, more pronounced minimum at shorter distance, in qualitative agreement with 
the RPA £ndingsi^>^>^ (for the top-fcc configuration this second minimum is not clearly 
visible). 

In Fig. 2, for the top-fcc case only, we also report the binding-energy curve obtained us¬ 
ing the DFT/vdW-WF2sl scheme where, however, at each graphene-Ni(lll) distance, the 
MLWFs used for the calculation of the vdW energy are just those obtained at the largest dis¬ 
tance (5.13 A) but rigidly shifted along the vertical 2; coordinate in such a way to be placed 
at the same level of the graphene layer (DFT/vdW-WF2sl-shift scheme). As can be seen, 
this curve is very different from that obtained by considering the MLWFs consistently gener¬ 
ated at each graphene-Ni(lll) distance; in particular it predicts that the equilibrium value 
(see Table II) is at a much larger distance, as expected for a conhguration dominated by 
standard long-range vdW interactions. This clearly demonstrates that the rearrangement 
of the MLWFs when graphene approaches the Ni(lll) substrate, as a consequence of the 
bond formation, is indeed crucial for a reasonable description of the energetics of the system. 
This observation is conhrmed by analyzing Fig. 3 where histograms are plotted showing the 
distribution of the WFCs along the z coordinate (perpendicular to the Ni(lll) surface), 
by also taking the relative spreads of the MLWFs into account as a Gaussian smearing. 
By considering, for the top-fcc conhguration, two different graphene-Ni(lll) distances, rep¬ 
resentative for the chemisorbed structure (corresponding to the binding energy minimum) 
and for the physisorbed structure (corresponding to the largest distance), respectively, the 
different distribution is evident. In fact, at large graphene-Ni(lll) distance there is a clear 
separation between the WFCs ’’belonging” to the graphene layer and those describing the 
electronic charge distribution of the Ni(lll) substrate (mainly localized at positions along 
2 ; corresponding to the locations of the metal layers), while, at the shorter equilibrium dis¬ 
tance, 3 WFCs, corresponding to MLWFs characterized by larger spreads (about 1.6 A) 
than those of typical MLWFs of free-standing graphene (about 0.8 A), are below graphene 
at an intermediate position between graphene and Ni(lll). By also considering (see Fig. 1) 
that 3 is just the number of C atoms located on top positions and thus able to form covalent 
bonds with the Ni substrate, within the simulation cell, this feature well illustrates the chem¬ 
ical bonding formation even without the need of a detailed analysis of the electron density 
distribution. Also note that the number of WFCs localized slightly above the topmost Ni 
layer increases from 2 to 6 as the graphene-Ni(lll) distance is decreased from 5.13 to 2.12 
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A, in line with the reportedi^ tendency of an electron transfer from the Ni(lll) substrate to 
adsorbed graphene. Moreover, the dashed curve obtained by the Gaussian smearing clearly 
suggests that a signihcant overlap between the MLFWs belonging to graphene and those 
of the Ni(lll) substrate occurs. Such an overlap is instead absent at the largest distance 
where evidently no chemical bond between graphene and Ni(lll) is formed. 


IV. CONCLUSIONS 

In conclusion, we have extended the applicability of our DFT/vdW-WF2sl method by 
investigating the interaction of Xe and graphene with the Ni(lll) transition-metal surface, 
the latter case representing a challenging scenario for a theoretical description. The results 
of our calculations, compared with available experimental data and other theoretical inves¬ 
tigations are rather good. We conhrm that the rare gas Xe atom is preferentially adsorbed 
on the top-site conhguration on the Ni(lll) surface too, similarly to what happens on other 
metallic surfaces. Moreover, our approach, based on the use of the MLWFs, allow us to well 
characterize the bonds between graphene and Ni(lll) which is clearly chemical in charac¬ 
ter is spite of the relatively small binding energy. This opens the way to applications of 
the method to investigate adsorption processes also on other interesting transition-metal 
surfaces such as Rh(lll). 
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TABLE I: Binding energy Ef, (in meV) and (in parenthesis) equilibrium distance R (in A) for 
Xe adsorbed on Ni(lll) in the top and hollow configurations, using different methods. In square 
parenthesis the adsorption energy E^ (in meV), which takes the Xe-Xe interaction into account 
(see text). Results are compared with available experimental reference data. 


method 

top 

hollow 

LDA 

-401 (2.51) 

-382 (2.17) 

PBE 

-19 (3.93) 

-18 (3.91) 

PBE+D 

-391 (2.76) 

-354 (3.25) 

rVVlO 

-284 (3.33) 

-271 (3.41) 

vdW-DF2 

-178 (3.80) 

-174 (3.84) 

DFT/vdW-WE2sl 

-283 (2.98) 

-229 (3.05) 

DFT/vdW-WF2sl 

[-369] 

[-316] 

expt.“ 

[-369] 

— 

expt.^ 

[-317] 

— 

expt.'^ 

[-191 ±9] 

— 














TABLE II: Binding energy (in meV per C atom, a positive value means that the configuration is 
unbound) and (in parenthesis) equilibrium distance R (in A) for graphene on Ni(lll), considering 
the top-fcc, bridge-top, and fcc-hcp configurations, using different methods, compared with available 


experimental and theoretical reference data. R is defined as the separation between the average z 
coordinate of the C atoms of the graphene layer and that of the topmost Ni atoms of the Ni(lll) 
substrate. 


method 

top-fcc 

bridge-top 

fcc-hcp 

LDA 

-196 (2.01) 

-210 (1.95) 

-35 (3.24) 

PBE 

+9 (2.17) 

-h9 (2.09) 

-5 (3.85) 

PBEsol 

-131 (2.08) 

-138 (1.99) 

-18 (3.49) 

DFT-D2“ 

-160 (2.09) 

-157 (2.05) 

-89 (3.02) 

DFT-D3^ 

-107 (3.25) 

— 

— 

rVVlO 

-81 (3.14) 

-80 (3.11) 

-75 (3.44) 

vdW-DF^^ 

-37 (3.50) 

— 

— 

vdW-DF2'^ 

-44 (3.68) 

— 

— 

009“^ 

-141 (2.07) 

— 

— 

M06L^ 

-44 (2.37) 

— 

— 

optBSGb-^ 

-112 (2.12) 

— 

— 

optBSSvdW® 

-71 (2.18) 

— 

— 

RPA^ 

-67 (2.17) 

— 

— 

RPA^ 

-70 (2.19) 

— 

— 

DFT/vdW-WF2sl 

-129 (2.12) 

-128 (1.99) 

-48 (3.64) 

DFT/vdW-WF2s 1-shift 

-29 (3.84) 

— 

— 

expt.* 

(2.11±0.07) 

— 

— 
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FIG. 1: Plane (upper panels) and side (lower panels) views of graphene on the Ni(lll) surface, in 
the top-fcc, bridge-top, and fcc-hcp configuration; the adopted simulation cell is shown in the plane 


views. 
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FIG. 2: Binding energy (per C atom) of graphene on Ni(lll) computed by the DFT/vdW-WF2sl 
method, in the top-fee (solid line), bridge-top (dashed line), and fee-hep (dot-dashed line) conhg- 
uration, as a function of the distance between the graphene layer and the Ni(lll) surface. The 
dotted line has been obtained using the DFT/vdW-WF2sl-shift scheme (see text). 
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FIG. 3: Histograms showing the distribution of the WFCs along the z coordinate (perpendicular 
to the Ni(lll) surface) for the top-fcc configuration for graphene on Ni(lll): (a) corresponds to 
the binding energy minimum at the 2.12 A distance, while (b) corresponds to the 5.13 A distance. 
The dashed curve has been obtained by a Gaussian smearing where the Gaussian widths are just 
given by the spreads of the MLWFs, while the vertical dotted lines denote the average position 
of the Ni atoms belonging to the different Niflllj lavers and (thicker line on the right) of the G 






































